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The influence of an external constant and homogeneous magnetic field H on the phase 
structure of the P-symmetric, chiral invariant 3-dimensional field theory model with two 
four-fermion interaction structures is considered. An arbitrary small (nonzero) magnetic 
field is shown to induce spontaneous violation of the initial symmetry (magnetic catalysis). 
Moreover, vacuum of the model at H ^ can be either P-symmetric or chiral invariant, 
depending on the values of the coupling constants. 



1 Introduction 

Magnetic catalysis is the phenomenon of dynamical symmetry breaking induced by external 
magnetic or chromomagnetic fields. This kind of influence of an external magnetic field H 
was for the first time observed in the study of the (2+l)-dimensional chiral invariant Gross- 
Neveu model |1|]. In that case, an arbitrary weak external magnetic field lead to dynamical 
chiral symmetry breaking (D^SB) even for an arbitrary small coupling constant [§, |3|. This 
phenomenon is explained by the effective reduction of space-time dimensions in the external 
magnetic field and, correspondingly, by strengthening the role of the infrared divergences in the 
vacuum reorganization Later, it has been shown that spontaneous chiral symmetry breaking 
can be induced by an external chromomagnetic field as well || ||, [7j . Moreover, basing upon the 
study of a number of field theories, it has been argued that magnetic catalysis of D^SB may 
have a universal, i.e. model independent, character ||. In a recent paper |J, a P symmetric 
3-dimensional Gross-Neveu model has been studied, and the external magnetic field proved 
to be a catalyst for spontaneous symmetry breaking again, this time for P parity breaking. 
Magnetic catalysis has already found its applications in cosmology and astrophysics [fLOR, and 



in constructing the theory of high-temperature superconductivity |], |TT[. It is evident that this 
phenomenon may be of crucial importance for various processes in the presence of an external 
magnetic field, where the dynamical symmetry breaking takes place, i.e. in the elementary- 
particle physics, condensed-matter physics, physics of neutron stars etc.Q 

As it is well known [[□[], the high temperature superconductivity is observed in antiferro- 
magnetic materials like La2Cu04, where conduction electrons are concentrated in planes formed 



by atoms of Cu and O. In recent experimental studies of this phenomenon [114] , it was found 
out that at temperatures much lower than the temperature of transition to the superconduct- 
ing state, the external magnetic field induces a parity-breaking phase transition. Thus, the 
problem of studying magnetic catalysis in high-temperature superconductivity becomes quite 



1 Possible applications of the magnetic catalysis were also discussed in recent publications |Tz 



urgent. The above arguments are first of all applicable to (2+l)-dimensional field-theoretical 
models of high-temperature superconductivity. 

In the present paper, we consider the influence of an external magnetic field on one of these 
3d models with a four-fermion interaction Lagrangian of the form 

N p N f~i N 

L = Y, ^ a + ^(E Maf + ^ a )\ (1) 

a=l a=l a=l 

where d = T^d^ , fields ip a are transformed according to the fundamental representation of 
the U(N) group introduced in order to employ the nonperturbative 1/iV-expansion method. 
Moreover, ip a is a four-component Dirac spinor for every value of a = 1,2, ...,N (correspond- 
ing indices are omitted). The action with the Lagrangian (HD is invariant under two discrete 
transformations, one of them being parity (P), and the other is chiral symmetry (these op- 
erations, as well as the matrix r are defined in the Appendix, where the algebra of V 1 matrices 
for the spinor reducible 4-dimensional representation of the Lorentz group is also presented). 
The model ([I]) at H = has already been studied in papers |]T3j| , and a corresponding gauge 



version of the model was considered in |T(|. This interest is mainly justified by the fact, that 



in a number of lattice models, describing dynamic effects in quantum antiferromagnetics in 
two spatial dimensions (e.g., high temperature superconductivity), going to continuum limit 
renders theories of the type ([!]). 

In contrast to simplest Gross-Neveu models, there are two different structures with four- 
fermion interaction in the Lagrangian ([I]). Hence, two different fermion condensates can exist 
in the theory: one of them, (ipip), is the order parameter for chiral symmetry, and the other 
one, (iprip), characterizes spontaneous breaking of parity. The condensate that provides for the 
lowest vacuum energy corresponds to the symmetry that is broken. The above arguments are 
fully applicable to the case of nonzero field H ^ 0. Therefore, the magnetic catalysis effects in 
the model (H) and in simplest Gross-Neveu models should be qualitatively different. 

Our main goal is to study characteristics of the magnetic catalysis in the P x r 5 -invariant 
model (P. At first, we are going to find out whether spontaneous symmetry breaking of the 
initial symmetry is induced by an external magnetic field. If this is the case, the question is 
whether both P and symmetries are broken simultaneously, or the vacuum of the model 
retains symmetry with respect to only one of the two discrete transformations. If the latter 
possibility is realized, then a question arises how the residual symmetry of the model depends 
on the values of coupling constants Gi^, etc. We demonstrate, that under certain conditions in 
the framework of the field theory ([]]), a parity breaking phase transition induced by an external 
magnetic field is possible. 



2 Phase structure of the model at H = 

Before considering the influence of the nonzero external magnetic field H on the vacuum of the 
model (0), its phase structure will be studied at H = 0. To this end, the auxiliary Lagrangian 

L = jidij) + a x {^) + (7 2 (VW0 + ^ + ^ (2) 

is introduced, where a\^, are auxiliary boson fields, and summation over indices a of the auxiliary 
group U(N) is implicitly performed here and in what follows. The field theories ([!]) and (@) are 



2 



equivalent, since by means of the equations of motion the fields 0^% can be excluded from (Q) 
and the Lagrangian ([1]) obtained. It is easily shown that the auxiliary fields are transformed 
with respect to discrete symmetries P and T 5 in the following way: 



P : a t (t,x,y) -xri(t, -x,y); a 2 (t,x,y) -> -a 2 (t, -a;,y); 
T 5 : o-i(t,x,y) -»• -at(t,x,y); a 2 (t,x,y) -»• a 2 (t,x,y), 



(3) 



i.e. cxi is a scalar field, and 02 is pseudoscalar one. Starting from the Lagrangian @, one can find 
the effective action of the theory, which has the following form in the one-loop approximation 
(i.e. in the leading order of the 1/N expansion): 



SefrK a 2 ) = -Nj d 3 x f ^- + ^-j - iSp In A, 



(4) 



where A = id + 01 + cr 2 r. Here, the fields <7 li2 are functions of space-time points. To obtain 
the effective potential of the model, the definition has to be used 



Vb (0-1,02) J d 3 x = -S e s(ai,a 2 ] 



<Tl,fT2=COIlSt 



(5) 



where boson fields are assumed to be independent of coordinates. With the help of (§) and (||), 
the following expression for the effective potential of the model ([!]) can be found in the leading 
order of the 1/N expansion: 



Vb (0-1,0-2) 



No 2 

2G7^" 2G 2 



No* 



2 r 

/ 

k=i j 



d 3 p 
(27rp 



ln(p 2 + M fe 2 ), 



(6) 



where M 12 = \<J 2 ± a±\, and integration is performed over Euclidean momentum. Integration 
in (pi) over the domain < p 2 < A 2 yields: 



V (a u a 2 )=Nj2 



k=l 



1 



2A 



( » 1 + 



G 



TV 



Ml 



(7) 



Now, in order to exclude the cutoff parameter A from (|?p, we introduce the renormalized 
coupling constants with the help of the following normalization conditions: 



1 d 2 V 
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N (da^ 2 






1 d 2 V 
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N (da 2 ) 2 


0-1=0 


= G~2 



2A 



7T 



2v x 

71 



2A 2u 2 



1 



7T 



(8) 



0-2=^2 



#2(^2) 

Relations (|8]) enable us to write the effective potential in terms of the ultraviolet finite quantities 



Vb (0-1,0-2) 



2 

k=l 



Ml 



— H 

2g k 6vr 



where (i = 1, 2): 



1 

ft 



ft(^i) 



2^ 

7T 



1 



1 

a 



(9) 



(10) 



and G c = n 2 / (2A). We remind that the coupling constants (^(z/j)) depend (do not depend) 
on the cutoff parameter A and do not depend (depend) on the normalization masses 2/j. By 
virtue of this, (0) leads us to the conclusion that the constants Oj depend neither on nor 
on A, i.e. the effective potential Vq is a renormalization invariant finite quantity (i.e., without 
ultraviolet divergences) . 

Thus, we have indeed demonstrated that the model (H) is renormalizable in the leading 
order of the 1/N expansion. The complete proof of its renormalizability is not intended in 
the present paper. However, we believe this to be true, relying on the results of the review 
paper |I7|, where the simplest 3-dimensional models with four-fermion interaction were proved 
to be renormalizable in the framework of the nonperturbative 1/N expansion method. 

It is well known that the phase structure of any theory is to a great extent determined by 
the symmetry of the global minimum of its effective potential. In order to avoid cumbersome 
formulas, we omit intermediate calculations, related to the study of the absolute minimum of 
the potential I'd (01,02) ® and present here the final result - the phase portrait of the model 
(|I|), depicted in Fig. 1 (similar analysis of the effective potentials for a whole series of the 
theories with the four-fermion interaction is carried out in || [Tj|, and, if required, it can be 
easily performed here). In this figure, the plane (01,02) of coupling constants fllOD is divided 
into five regions, in each of which one of the possible phases of the model is realized: A, B, 
C. In the phase A, where g>x,2 > 0, corresponding to small values of the coupling constants 
G\,2 < G c , the global minimum of the function Vo(0i,02) is at the origin (o~i = 02 = 0), and 
hence, the ground state of the system in this phase is P and Ts symmetric. In phases B and C, 
where at least one of the coupling constants Gi is greater than G c , the point of the potential 
global minimum is correspondingly at o\ = 0, cr 2 = — 7r/g 2 and o\ = —n/gi, a 2 = 0. Hence, 
with the aid of (^), one can see that the vacuum is Ts invariant in the phase B, and the parity 
is spontaneously broken. The ground state of the phase C of the model is P invariant, however 
the chiral symmetry is spontaneously broken in this case. 



3 Effective potential at H ^ 

Before starting to investigate the vacuum of the 3-dimensional model ([]]) in the external constant 
magnetic field, let us work out an expression for the effective potential at H ^ 0. The Lagrangian 
of the model in terms of auxiliary scalar fields 01,2 (x) takes the form 

L H = $(id + ei)V + 01 {jH>) + 02OMO + ^ + (II) 

where A = A^V 1 , e is the (positive) charge of fermions, and the vector-potential of the constant 
external magnetic field H has the form Ao,i = 0, A2 = X\H. Using the standard technique of 
the 1/N expansion |TJ, one obtains, in the leading 1/N order, the following expression for the 
effective potential: 

Ir , , A^0? A^0| iN^ 

Kfr(0i,0 2 ) = ^ + ^ + — Spin A H . (12) 

Let us note, that here A# = id + eA + o\ + r<7 2 is an operator that acts only in the spinor and 
coordinate spaces, v = J d 3 x, and fields 0j do not depend on space-time points. 

First, let us turn to the causal Green's function of the operator Ah- It can be presented in 
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the following form: 



{n} 

- t) £ ^(f,^U/9(^0- (13) 

M 

Here {n} = (z, n, fc), where i — 1, 2; n = 0, 1, 2, ...; A; is a real number — oo < k < oo. Moreover, 
■0/^1 are positive- and negative-frequency orthogonal normalized solutions of the Dirac equation 
Aff^ = 0, which in the case a± + a 2 > 0, o 2 < o~i have the form (T is a symbol of transposition): 

^inl = exp(=pe n t + ifcc 2 ) X 



x I W — i — fcn-M^O.O.O 



(14) 



x (o,o,y — ^(x,),^ ^ K_ 1>k ( Xl ) 



where 



(etf) 1 / 4 
(2 n n\ 



Kk(xi) = /n l^ y/2 ex P (-e/2)Hn(0, (15) 



e n = y (ci + + 2eHn, e n = J (cri — a 2 ) 2 + 2eHn } H n (£) are the Hermite polynomials, 
£ = \/ eH [xi — k/{eH)\ . Functions (|15|) satisfy conditions 

J d xi K,ki x i) = J dk h l,k(xi) = 1- (16) 



Moreover, in formulas (13), it is assumed that h_i^(xi) = 0. 

Now we have obtained everything necessary for calculating the value A(a 1 ,a 2 ) = 
zf _1 SplnA^. It is evident, that 



dA(a 1 ,a 2 

— IV 



d<7\ 



w _1 SpA^ = iv' 1 f d 3 x ^](A H 1 ) aQ (x,t;f,t). (17) 

^ a=l 



Taking formulas (|T~3D — (|16D into account, after rather tedious calculations, one obtains from this 
expression (similar calculations were performed in the case of the simplest Gross-Neveu model 
in §): 

dA(a 1 , cr 2 ) {(Ti + cr 2 )eH °f ds 



/ds 
^exp(-^ 1 + a 2 ))cth( e if S ) 

v 

(a 1 — o 2 )eH f ds . . s9s , . TT . 

/— exp -s d! - a 2 2 eth ers . 
J \ s 



v 

Integrating both sides of this equality over o\ within the limits o\ and 00, one obtains 

eH 2 7 ds 

A(a u a 2 ) - A{oo, a 2 ) = — ^ £ J — 2 exp(- S M 2 ) cth(eHs), (19) 
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where expressions for M, are presented in (y). Note, that formula flT9|) was obtained with the 
assumption that a 2 > a 2 < &\- However, one can directly prove that ([19]) is also valid for 
arbitrary relation between o\ and a 2 . 

The value of dA(a\, a 2 )/da 2 can be found in a similar way. One can integrate it over a 2 
within the limits a 2 and 00, and then the expression obtained can be compared with (|19|). As 
a result, we obtain the equality A(oo, a 2 ) = A(ai, 00), which leads to the conclusion that, due 
to mutual independence of variables <i\ and a 2 , the value of A(oo,a 2 ) does not depend on a 2 , 
i.e. it is constant. 

Taking this into account, due to equation fllTf), one obtains from (|T2| ) the following expression 
for the effective potential (up to an inessential constant): 



,r , \ Na? Nal 2 e H r ds , , ~. , . rT , 
VWl,*) = ^ + ^ + N f^2 ^exp(- S A4 2 )<=th(eff S ). 



(20) 



This expression evidently contains UV divergent integrals. Upon some transformations in (|20|) , 
this divergence can be localized in the part of the effective potential that corresponds to H = 
contribution Vo(o"i,cr2): 



2 eH 7 ds 
Kfr(ai,a 2 ) = V (a 1 ,a 2 ) + N^2 — 

^3/2 J ~3/2 eX P( — 
k=l n 



cth(eHs) 



eHs 



(21) 



Here 



Vo(ai,a 2 ) 



Na\ Na 2 2 A 2 

H + N^2 

2G\ 2G 2 k=1 



8vr 3 /2 



ds 
7572 



exp(-sM 2 ; 



:22i 



It is easily seen that only the first term in (|2"T|), defined in (|2"2"D, is singular. Up to a certain 
additive constant, it is equal to potential (^). Upon renormalization, it evidently takes the form 
of (||). The following expression for V e s, resulting from ( |2"TD after integration over s (see 
will also be used: 



V eS (a u a 2 ) = Nj2 



k=l L 



a\ | eHM k (2eH) 
2g k 47r 47r 



3/2 



-c 



1 Ml 



2' 2eH 



(23) 



where C( s i x ) 1S the generalized Riemann Zeta-function [20]. 



4 The phenomenon of magnetic catalysis 

Now, let us consider spontaneous symmetry breaking in the initial model ([!]). To this end, 
properties of the global minimum of the potential (^) under symmetry transformation (|J) 
have to be investigated. In this case, positive values of parameters g^ 2 , i.e., correspondingly, 
the region of small values of the bare coupling constants Gi t2 , where the P x r 5 -symmetric 
phase of the model exists at H = 0, are of special interest. The question is what kind of 
symmetry the vacuum possesses at H 7^ in the case under consideration, i.e. at gi > 0, 
g 2 > 0. To find an answer to this question, it is useful to consider a simplified situation. 

Special case of model ([I]). First, we illustrate the magnetic catalysis phenomenon by 
considering the example of a simplest model that follows from (|l|) when G\ = G, G 2 = 0. 
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The corresponding effective potential Vh 9 (<j) can be obtained from (23), assuming o-i = 0, and 
introducing the notations <7i = a, g± = g: 

_ _ Na 2 NeH\a\ N{2eHf' 2 ( 1 a 2 \ 

It is evident that function Vff 9 (cr) is symmetric under the transformation a — > —a. Hence, in 
order to find its global minimum, it is sufficient to study V/j s (cr) only in the domain a G [0, 00). 
Here, the stationary equation has the form 



dV Hg (a) _ Na NeH NaV2eH (l 



a 



2 



+ ^ S =°- ( 25 ) 



cV g 2ti 2tt s ^2' 2e# 

This can be obtained with account for the fact that 

d((s, x)/dx = —s((s + l,x). 
When cr — > 0, the following expansion for the ^-function is known 



C (i = + const + od^l/v/M). (26) 



Substitution of (p6|) into a LHS of (p5|) gives 

dV Hg (a 



da 



(27) 



2tt 



i.e. the point a = is not a solution of the stationary equation (E3|). Moreover, since Vff g (cr) = 
Vh 9 (— 0"), two more consequences follow from (27): 1. The point a = is a local maximum of 
potential (^4]). 2. The first derivative of function Vu g {cr) does not exist at the point a = 0. (At 
the same time, the effective potential for H = is a differentiate function in the whole axis 
a = 0.) Consequence n.l and the fact that \im\ (T \_ >00 VHg{o') = +00, imply that there exists a 
point ao(H,g), where the potential Vjig(cr) has a nontrivial, i.e. nonvanishing global minimum. 
This means that chiral symmetry of the simplest model in question is necessarily spontaneously 
broken in the presence of an arbitrary weak external magnetic field, when the coupling constant 
G is arbitrarily small (constant g is positive). Thus, dynamical breaking of chiral invariance, 
induced by an external magnetic field, takes place (magnetic catalysis). More detailed analysis 
of the phase structure of the simplest Gross-Neveu model in an external magnetic field can 
be found in || ||, || [21], where the influence of temperature, chemical potential and nonzero 
space-time curvature on the magnetic catalysis is also studied. [] 

In what follows, more detailed information on the properties of the effective potential (p4|) 
will be needed. Assume that g > and eHg 2 1, i.e. the magnetic field is weak. In this case, 
with the use of expansion ( p6|) in equation fl25|) , one can easily obtain the following asymptotic 
behavior of the global minimum point a (H, g) of potential Vh 9 (cf) at eH — > 0, being a solution 
of the stationary equation (|25|) : 

a (H,g) = egH/(27r) + .... (28) 



influence of such factors as temperature, chemical potential, nontrivial metrics and topology of the space 
(disregarding magnetic field) on the ground state structure in the simplest (2+l)-dimensional Gross-Neveu 
model is considered in [E2I. 
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We now consider the quantity VH g (ao(H, g)). It is evident that for arbitrary values of H 
and g the following relation holds: 



dV Hg (a (H,g)) _ da (H,g)dV Hg (a) 



dg 



dg 



da 



+ 



dV Hg {a) 



o=o Q {H )g ) 



dg 



(29) 



a=<r {H, g ) 



The first term in the right-hand side of the equation vanishes in virtue of the fact that ao(H, g) 
satisfies the stationary equation (|25|). Thus, we have 



dV Hg (a {H,g)) _ dV Hg (a) 



dg 



$9 



2g 2 



< 0. 



(30) 



It follows from (|30|) that VH g {ao(H, g)) is a monotonically decreasing function of parameter g. 

Magnetic catalysis in the general case. We now consider spontaneous symmetry brea- 
king in model ([!]) under the influence of an external magnetic field in the general case, when two 
coupling constants G\ i2 are arbitrary First of all, we will demonstrate that the vacuum of the 
model at H ^ is not P xT 5 symmetric any more, i.e. the point a\ )2 = does not correspond 
to the global minimum of the function V e s(a\, a 2 ). To this end, we study first derivatives of 
this function. They have the following form in the domain (a\ + a 2 ) > 0, a k > a\ (here k ^ I): 



1 dV eS (ai,a 2 ) 



N 



da k 



9k 



eH 
27 



*i 



a 2 )V2^H (I (at + a 2 ) 
4vr C \2' 2eH 



{a k -a l )V2eH fl (a 1 - a 2 ) 2 " 
4vr H 2 ' 2eH 



1 dV eS (ai,a 2 ) 



N 



da. 



ai_ _ (a 1 + a 2 )V2^H ( 1 {ax + a 2 ) T 
gi 4tt U' 2eH 



+ 



(a k -a l )V2eH ( I 
+ 4tt M 2' 



Pi 



0"2 



2eH 



(31) 



If we put k — 2, I — 1, o\ — first, and then k — 1, I 
easily see, with the help of expansion (|26l), that 



2, a 2 = in these formulas, one can 



dV P 



eff 



da 2 



%2 ' 



NeH 
2tt 



< 0; 



eff 



dai 



21 

(T 2 : 



2tt 



< 0. 



(32) 



Relations (|32|) imply that, when H ^ 0, effective potential (|23| ) decreases at the point a\^ 2 = 
in the positive directions of axes a\ and a 2 . Hence, this point can not correspond to a global 
minimum of the potential V e s (a\, a 2 ), and an arbitrary weak external magnetic field induces 
spontaneous breaking of the initial PxT 5 symmetry for arbitrary coupling constants G\ i2 . 

Nondifferentiability of functions (p3|) on the line a\ = a 2 is another important consequence 
of formulas (31). This results from the fact that each of the partial derivatives (31) takes 
different values, when approaching this line from above and from below. Moreover, since the 
function (^) is symmetric against each of the transformations ai — > — o^, a 2 — > —a 2 , it is 
evidently nondifferentiable in the plane of variables (ai, a 2 ) on the lines o~i = a 2 and a\ = —a 2 . 
We note that at H = the effective potential is different iable in the whole plane (a±, a 2 ). The 
function (^) is even in each of its variables and hence, we will further study its global minimum 
only in the domain o\ > 0, a 2 > 0. 



S 



Consider the influence of an external magnetic field on the symmetric phase of the model in 
more detail, i.e. we assume that g\ and g 2 are positive. It follows from (|I0D that bare coupling 
constants are sufficiently small in this case, i.e. Gi < G c = 7r 2 /(2A). As it was shown above, 
the magnetic field destroys the P x T 5 symmetry of the ground state of the model, i.e. magnetic 
field is a catalyst of the spontaneous breaking of this symmetry. In order to study magnetic 
catalysis in more detail, we first assume that the magnetic field is so small that eHgf <C 1, 
eHg\ <C 1. Then only two points satisfy the system of stationary equation dV^jda^ = (where 
k = 1,2): i. ( 0io,O) and ii. (0, 020), where a i0 = a (H,gi), and a (H,gi) are the solutions of 
equation (|25|), where g = gi- It is evident that the global minimum of potential V e s (01, 02) can 
be situated only in one of these two points on the (cj., 02)-plane, where the following relations 
hold 

Kff (010,0) = V Hgi (a (H, gi )); K ff (O,0 2o ) = V Hg2 {a (H, g 2 )). (33) 

In the right hand sides of these equations, the minimum value of function (|24]) at g — g\ and 
9 = 92 stands. According to ( ^9|) and (p0|), the values of the function in the right hand sides of 
these equations decrease monotonically with parameter g. Therefore, at g\ < g 2 (g 2 < gi) the 
effective potential fl2l"D is smaller at the point ii (i) than at the point i (ii). 

The final conclusion about the position of the global minimum of potential (^) can be 
drawn only after studying this function on the line L = {(01, 02) : u 2 = 01}. On this line, as it 
was shown above, the potential is not differentiable, and hence, there exists a probability that 
the global minimum of the potential ^(01,02) is situated just here. Consider in more detail 
such a possibility, and to this end, we study the function ( p3l) on the straight line L where it 
has the form 

K cff (0, 0) = V[a) = — + — - — . (34) 

In this expression, g = g\g 2 j (gi + g 2 ), and the function VH g (cr) is defined in (0). The stationary 
equation for ^(0) has a unique solution O , which is the absolute minimum of this function 
(this is demonstrated in the same way as for the potential (|24|)). It can be shown that the 
quantity 00 is also small in the regions where the magnetic field is weak, and moreover 

00 = a (H,g) + o(eH) = eH~g/(2n) + o(eH). (35) 

Therefore, to estimate the magnitude of ^(00), one can use the relation 

F(x) =F(0) + F'(0 + ) -x + o(x), at x 0+ , (36) 

where F(x) is an arbitrary function differentiable at the right from the origin, i.e. there exists 
a finite limit F'(0 + ) = \im x ^ 0+ (F(x) —F(0))/x. Substituting a = a in (|34]) and using equality 
(|36|), one can easily obtain the following value of the function V(cr) at the point of the absolute 
minimum at H — > 0: 

T~r,~ N t r fn , dV H g(0) _ A^0Q 

« «-« (37) 

The second relation in this expression can be obtained with regard to formulas (B^) and (|35|). 
In the same way, one can estimate the quantities in formula 



V Hgi (a (H, 9i )) « V Hg M - ^J 9 * - (3? 
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Comparing expressions fl3?D and fl38|) , one can easily conclude that, in the set > 0, the 
potential ( pT|) can not reach its minimum on the straight line L. Therefore, for g 2 > g± the 
global minimum of the potential is at the point (0,cr 2 o), and for g x > g 2 it is at the point 

K),o). 

Thus, when an arbitrarily weak external magnetic field is applied at the P x r§ symmetric 
A-phase of the model (for which two parameters g\^ are positive), ground state of this phase 
is destroyed. For g 2 > g\, the vacuum maintains T 5 invariance, and phase A transforms 
to B phase of the theory. At g\ > g 2 the ground state is only P symmetric, and phase C 
spontaneously appears instead of phase A. As far as o- (H,g) — > at H — ► 0, in both cases, 
magnetic field induces the continuous second order phase transitions from phase A to phases 
B or C correspondingly. 

Unfortunately, for the magnetic field such that eHgf ~ 1 or eHgf ^> 1, we were unable to 
prove this fact by means of the analytical technique. However, numerical calculations performed 
for a rather wide interval of values of H confirm the above conclusion. For example, Fig. 2 
depicts the function V e s(ai, <r 2 ) at H = 40e 3 , g± = e~ 2 , g 2 = 20e~ 2 (here, evidently eHgf ^> 1). 
In this case, the global minimum of the potential a± = 0, a 2 = 2.77 e 2 corresponds to chiral 
symmetric phase of the model (parity is spontaneously broken). 

It should also be remarked that an external magnetic field acting upon phases B and C (see 
Fig. 1) makes their ground states more stable. The fact is that there exists a rather strong 
fermion-antifermion attraction even at H = (at least one of the bare coupling constants Gi is 
greater than G c ), leading to formation of one of the condensates ipij} or iprip, i.e. to spontaneous 
breaking of the initial symmetry. An external magnetic field acts in favor of further increasing 
the condensate values, i.e. it stabilizes these phases. Since in these cases the magnetic field is 
not the cause (catalyzer) of the spontaneous symmetry breaking (a superstrong attraction of 
fermions and antifermions being the actual cause), details of these considerations are omitted. 
For completeness, however, the phase portrait of the model at arbitrary nonzero values of the 
external magnetic field and coupling constants g\ )2 is depicted in Fig. 3. 

5 Conclusions 

In the present paper, the influence of an external magnetic field on the 3-dimensional P x 
symmetric model ([!]) is considered in the case, when bare coupling constants Gi ;2 < G c take 
sufficiently small values (positive values of constants gi j2 ) . The initial symmetry is unbroken, 
if H = 0. Arbitrary small external magnetic field was shown to induce spontaneous breaking 
of P x T 5 symmetry (phenomenon of magnetic catalysis). Moreover it turned out, that P 
symmetry of the model is broken, while chiral symmetry remains intact, if g 2 > gi (this is 
phase B of the theory). When g\ > g 2 , the external magnetic field induces the second order 
phase transition from phase A to phase C, whose vacuum is P symmetric, though having no 
chiral invariance. 
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Appendix 

Algebra of T-matrices for the 3-dimensional Lorentz group 

Two-component Dirac spinors realize an irreducible 2- dimensional representation of the Lorentz 
group transformations of the 3-dimensional space-time. In this case, the 2x2 matrices 7 have 
the form 

7° = -s=(j _° 1 ),y = ia 1 =(° * ),-y 2 = ^=(_? 1 j)- (A.l) 
These matrices have the properties 

Tr(YY) = 2gT\ [Y, Y] = -2ie^ a ; YY = -ie^la + <T , (A.2) 
where g^ v = = diag(l, —1, —1), 7 Q = g a pY 1 £ 012 = 1- Moreover, 

Tr{YYY) = -2ie^ a . (A.3) 

Recently, one frequently exploits a reducible 4-dimensional representation of the Lorentz 
group, and the spinors i[} k from (1) transform according to this representation: 

*-(£)■ (A - 4) 

In (A. 4) ipi,ip2 are two-component spinors. The gamma matrices that correspond to this 
representation have the form T M = diag(7^, — 7^), where 7^ are given in (A.l). It is easily 
demonstrated that (fi, v = 0, 1, 2): 

Tr(rT") = Ag^- m u = o» v + g^; 

a v» = I[r", Y v \ = diag(-i£^ 7a , -ie^Ja). (A.5) 

The dimensionality of the matrix algebra that acts in the four-dimensional spinor space is equal 
to 16, and its generators are (/i = 0, 1, 2) and the matrix T 3 (anticommuting with them): 



I 

1 



(A.6) 



where J is a unit 2x2 matrix. There exists one more matrix T 5 , which anticommutes with all 
matrices r M and T 3 : 

r^(?-'). ,a,) 

We note that matrix r entering ([!]) has the form 

r-(j^). (A.7, 

In the space of four-dimensional spinors (A. 4), one can introduce the parity (P) and chiral r 5 
transformations, those in terms of two-component spinors ip\^ look like 

P : i/j lk (t,x,y) — > Yip2k(t } -x,y) } ifa k (t,x,y) — ► Y^ik(t, -x, y), 

T 5 : if)i k (t,x,y) — > -iif) 2k (t,x,y), if) 2k (t,x,y) — >• ii[>i k (t, x, y). (A.8) 

Discrete transformations (A.8) can be easily rewritten with the use of the four- component 
spinors (A. 4) 

P: ^{t,x,y)^iT l Y^{t,-x,y)- r B : V - Tfy (A.9) 
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Figure 1: Phase portrait of the model in terms of gi and g 2 &t H = 0. Here / = {(gi, g 2 ) '■ gi = 

92}- 




Figure 2: Diagram of the function (<5"i, 02) = V e g((Ti, o"2)/e 6 at H = 40e 3 , g\ — e 2 , g>2 = 
20e~ 2 . Here <7i = ai/e 2 , <t 2 = o" 2 /e 2 . 
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Figure 3: Phase portrait of the model in terms of g\ and gi at H ^ 0. Line / defined on Fig. 1. 



